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Abstract. Disordered quantum networks, as those describing light-harvesting 
complexes, are often characterized by the presence of antenna structures where the light 
is captured and inner structures (reaction centers) where the excitation is transferred. 
Antennae often display distinguished coherent features: their eigenstates can be 
separated, with respect to the transfer of excitation, in the two classes of superradiant 
and subradiant states. Both are important to optimize transfer efficiency. In absence 
of disorder superradiant states have an enhanced coupling strength to the RC, while 
subradiant ones are basically decoupled from it. Disorder induces a coupling between 
subradiant and superradiant states, thus creating an indirect coupling to the RC. 
We consider the problem of finding the maximal excitation transfer efficiency as a 
function of the RC energy and the disorder strength, first in a paradigmatic three- 
level system and then in a realistic model for the light-harvesting complex of purple 
bacteria. Specifically, we focus on the case in which the excitation is initially on a 
subradiant state, showing that the optimal disorder is of the order of the superradiant 
coupling. We also determine the optimal detuning between the initial state and the 
RC energy. We show that the efficiency remains high around the optimal detuning 
in a large energy window, proportional to the superradiant coupling. This allows for 
the simultaneous optimization of excitation transfer from several initial states with 
different optimal detuning. 

PACS numbers: 71.35.-y, 05.60.Gg 

Keywords: quantum transport in disordered systems; open quantum systems; light¬ 
harvesting complexes. 
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1. Introduction 

Photosynthetic bacteria utilize antenna complexes to capture photons and convert the 
energy of the short-lived electronic excitation in a more stable form, such as chemical 
bonds. After absorption, the energy is transferred to a complex, called reaction center 
(RC), where it initiates electron transfer, resulting in a membrane potential. This very 
efficient transfer occurs on a time-scale of few hundreds of picoseconds and on a length- 
scale of few nanometers, so that coherent quantum dynamics can enter the play, as recent 
experiments seem to prove [1]. Quantum coherence can enhance transport efficiency 
inducing Supertransfer and Superradiance in light-harvesting complexes [21 El II]- On 
the other hand, quantum coherence can also be detrimental to transport, as Anderson 
localization and the presence of trapping-free subspaces [B] show. 

Superradiance [3 El El, as viewed in the context of both optical fluorescence 
unmimg and quantum transport in open systems na m ns EH ES], is not solely a 
many-body effect. Single-excitation superradiance is a prominent example of genuinely 
quantum cooperative effect nzi, relevant in natural complexes, which operate in the 
single excitation regime since solar light is very dilute. 

Natural complexes are subject to a noisy environment with different correlation 
time-scales (if compared to the excitonic transport time): (i) short-time correlations, 
giving rise to dephasing (homogeneous broadening, as considered e.g. in [T8l [TQ] 1201121] ) 
and (ii) long-time correlations, producing on-site static disorder (inhomogeneous 
broadening, as considered e.g. in [221125] ). The role of environmental noise is twofold: 
on one hand, it can help transport since it destroys the detrimental coherent effects, 
leading to noise-enhanced energy transfer, i.e. the existence of a maximal efficiency 
at some intermediate noise strength, as found in the last decade by various groups 
[2n [25l fTO] 120] [6l [261 EHl E3 Ell [291 [211130] . On the other hand, it can suppress the 
beneficial coherences leading to a quenching of Supertransfer na. It is thus essential to 
consider this non-trivial interplay when dealing with realistic light-harvesting complexes, 
as it will be shown below. 

The typical structure of bacterial photosynthetic complexes displays a RC placed at 
the center of the light-harvesting complex I (LHI), with the chromophores arranged on a 
ring and surrounded by other ring structures (called LHII) acting as peripheral antennae. 
The LHI-RC structure describing light-harvesting complexes produces a distinguished 
feature: the eigenstates of the peripheral ring structure can be separated into two classes, 
superradiant and subradiant states, with respect to the transfer of excitation towards the 
RC. In absence of disorder, the superradiant states are coupled to the RC with a coupling 
amplitude proportional to \/iV, where N is the number of chromophores in the ring, 
while subradiant states are basically decoupled from the RC. In presence of disorder, 
the subradiant states can be coupled to the superradiant ones and, as a consequence, 
only indirectly coupled to the RC states. Note that the subradiant subspace, previously 
studied by the authors [TB], has been also analyzed in the literature under the name of 
trapping-free subspace [6]. 
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Figure 1. In panel (a) the trimer model is shown in the site basis: two sites connected 
with the RC with equal coupling Hj-c and between them with coupling Vi. In panel (b) 
the trimer model is shown in the superradiant-subradiant-RC basis, see text. 


Here we want to discuss the optimization of excitation transfer efficiency from 
peripheral states in networks displaying the Ring-RC structure, with particular reference 
to subradiant initial states. Indeed, since disorder is only detrimental for the 
superradiant state, we expect that the optimal condition will be mainly determined 
by the behavior of subradiant states. Nevertheless, the role of superradiant states is 
very important [15] and it will be the subject of future investigations. 

Even if in natural systems the role played by dephasing is crucial, here we will focus 
only on the interplay of static disorder and energy landscape (energy detuning between 
the initial state and the RC) on transport efficiency. 

To pursue our goal, we consider in section a paradigmatic model: a trimer 
displaying a superradiant and a subradiant state, on which the excitation is initialized, 
and an acceptor state (the RC), where the excitation can be trapped, see hgure 
In absence of disorder, while the subradiant state is decoupled from the RC, the 
superradiant state has an enhanced coupling Ij-c to the RC. We discuss the maximization 
of the average (over disorder) transfer efficiency with respect to both the detuning hErc 
between initial subradiant state and RC and the strength W of disorder. We have found 
two different regimes depending on Id-c and the energy gap A between the subradiant and 
superradiant states. While the optimal noise is always of the order of I^c, the optimal 
detuning occurs at resonance with the initial energy, (jErc = 0, only for Id-c > A. Another 
optimal detuning, 5E.rc = should be used in the case Id-c ^ A. 

In section we show the relevance of our Endings in a more realistic model by 
analyzing the LHI-RC complex of purple bacteria as presented in EH. After a reduction 
of the whole system to simpler systems capable to describe its main features, we prove 
that the optimal conditions for the excitation transfer from a room-temperature Gibbs 
distribution of ring states can be estimated by applying our analysis. The optimal 
disorder strength we find falls within the estimated range of natural disorder strength. 
We conclude the paper (section with a summary of the main results and of future 
research directions. 
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Many transport properties of systems displaying the Ring-RC strnctnre can be 
nnderstood from the analysis of the simplest ring, i.e. two conpled sites, connected 
with a central site representing the RC, as depicted in £gure[^. The Hamiltonian (in 
site-basis) for snch trimer model can be written in matrix form as follows: 

/i?i —irfl/2 H flj-c \ 

I H E2—iTQ^/2 fli-c j (1) 

\ f^rc Ej-c — iiv-\- k)/2 ) 

where the action of the environment (static disorder) has been taken into acconnt by 
choosing the energy levels Ek, k = 1,2, as Ganssian random nnmbers with mean zero 
and variance W"^, and no disorder has been added to the RC site. 

The loss of excitation through the RC has been described by the non-Hermitian 
term —m/2. Throughout the whole paper, k is assumed to be small with respect to the 
other coupling parameters. This choice is consistent with the realistic photosynthetic 
model which will be considered in the next section. In order to make a close comparison 
with realistic systems and following a standard procedure [1 El OEl H [IS], we also 
introduced the diagonal non-Hermitian terms —iTf[/2, with fluorescence constant Ta 
much smaller than any other energy scale (C, Crc, kk, k), which represent the loss of 
excitation from each site due to recombination. 

It is convenient to move from the site-basis to the subradiant-superradiant-RC basis, 
by dehning the states 


\SUB) = ^{\l)-\2)) , 


\SR) = ^{\l) + \2)) , 


from which the new Hamiltonian El easily follows, see also hgure [^: 

f—VL — iT^/2 X 0 \ 


H = 


X 


C —irfl/2 


^/2VLr 


( 2 ) 


\ 0 V^Crc y - i{Tf[ + k )/2 / 

where the two Gaussian random variables x = {Ei — E2)/2 and y = E^.^ — {Ei -|- E2)/2 
have been introduced and are such that 


(a;)=0, {y)=E,,, {x^) = {y^) - El = W^/2 . 

In this basis, the coupling between the subradiant state and the RC vanishes, whereas 
the coupling between the superradiant state and the RC is enhanced and it is given 
by the matrix element l^c = V^Ci-c. Note that in a ring of N sites, the superradiant 
coupling Vrc typically scales as a/ZV. The structure of the Hamiltonian, depicted in 
hgure [^, implies that the excitation transfer from the subradiant state to the RC can 
only be mediated by the superradiant state, through the random coupling x. 
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Since we are interested in modeling light-harvesting complexes, a most relevant 
quantity is the efficiency, at the time t, of energy transfer from the system into the RC. 
Given an initial state |ft^in), it is dehned as 

I (RC I , (3) 

\ Jo / w 

and it represents the probability of escaping out of the system up to the time t. In the 
above dehnition the brackets {.. .)w indicate the average over disorder. 

In numerical simulations we always consider the efficiency at a time t > h/Tfi. Note 
that the efficiency strongly depends on the time t at which it is computed for t < h/Tfi, 
while it reaches a stable asymptotic value 1]°° for t > h/Tfi, thus motivating our choice. 
As Ffl —)■ 0, the asymptotic value of the efficiency is ri°° = 1 for any choice of parameters. 
Hereafter we will measure energies in cm“^ and times in ps, that corresponds to setting 
l//i ~ O.Obvr cm/ps. 

It is clear from ([^ that the energy transfer efficiency is strongly dependent upon 
the initial state, a feature also studied in |23]. Indeed, if we start from the superradiant 
state we are in a situation in which there is an enhanced direct coupling to the RC at zero 
disorder. Thus, one might think that the best situation occurs when the excitation is 
on the superradiant state set at resonance with the energy of the RC. In this situation 
disorder is only detrimental to transport, since it tends to destroy the superradiant 
coupling [TOl [Is] and it moves the system out of resonance. On the other side, the 
excitation in natural complexes is usually spread also on subradiant states, due to the 
presence of a thermal bath. Since a subradiant state is not directly coupled to the RC, 
it is only through the action of disorder that the excitation can be transferred from the 
initial state to the RC (i.e. for IR = 0, we have rj = 0). 

We will focus our attention on this non-trivial case (see hgure &) , in which an 
initial excitation is on the subradiant state at energy Ei^i = —0, coupled via x to the 
superradiant state at energy Esr = = i?in + A, which is further coupled to the RC 

with tunnelling amplitude Re = V^flrc- Our aim is to hnd the system conhguration 
which maximizes the average transfer efficiency at f > ^/Tfj. Fixing the energy gap 
A = 2f2 between the superradiant and the subradiant states and the superradiant-RC 
coupling Rc, and assuming n and Fa to be perturbative quantities, we are left with two 
independent parameters to be tuned to achieve the maximal efficiency: the subradiant- 
RC detuning 6E,.c = E^c ~ Em and the strength W of the random coupling x. 

To pursue our goal we will hrst analyze a fully deterministic model obtained 
replacing the stochastic terms x and y in equation ([^ with deterministic parameters as 
follows: 

/ -E-in — *rfl/2 X 0 \ 


^^det 


A E[^ -|- A — iFfl/2 Rc 


(4) 


Rc 


\ 0 


Ere ~ *(rfl -|- k) j2 / 
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Figure 2. Average transfer efficiency g at fluorescence time tfl = S/Ffl, computed 
starting from the subradiant state \SUB), plotted in the left panel as a function of 
the subradiant-RC detuning 5E^c and of the deterministic subradiant-superradiant 
coupling X, and in the right panel as a function of JA'rc and of the rescaled disorder 
strength W = Wl\/2. The blue cross Indicates the estimate ([^ for the optimal transfer 
conditions. The dashed curves indicate the resonances determined in (|^. The value of 
the parameters are A = 0 cm“^, Vrc = 1 cm“^, k = 0.01 cm“^, and Ta = 10“'* cm“*. 


In particular, we will discuss which pair , Xopt) of values of the subradiant-RC 

detuning and coupling strength produce the maximal efficiency. Subsequently, we will 
investigate how disorder affects those results. 

In what follows, we will hrst analyze the results of the deterministic model and then 
the results in presence of disorder. In hguresj^ and the results of the deterministic 
model are shown in the left panels, while those in presence of disorder, are in the right 
panels. 

2.1. Analysis of the deterministic model 

According to our previous assumptions, the behavior of the efficiency p in the {X, dEj-f)- 
plane depends on the ratio between the two system parameters A and Kc- Since it is 
relevant for natural systems, we hrst consider the case A = 0 , which corresponds to two 
uncoupled ring sites in the trimer model (f2 = 0) equally connected to the RC. In this 
case, the sole energy scale of the system is I^c- We will subsequently consider the effects 
of a hnite gap A 7 ^ 0. 

2.1.1. The zero-gap case With A = 0 our model corresponds to a tight-binding chain 
of three sites, the hrst (subradiant) is coupled via X to the second (superradiant), with 
equal energy, which is coupled to the RC with strength Id-c R can be easily checked that, 
if we initially excite the subradiant state, the probability Prc(^) of hnding the excitation 
on the RC (which is a periodic function of time if we neglect the non-Hermitian terms 
in (|^) can reach the maximal value of 1 in the shortest time when X = Vrc and 5Ere = 0 , 
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(a) 


(b) 


V 


X 


-E.v 



Figure 3. Different schematic representations of the energy levels in the deterministic 
trimer model. (a) The subradiant-superradiant subsystem (framed), where the 
coupling X produces the dressed levels with energies e±, is coupled through Vi-c 
to the RC state at energy i?rc- (b) The initial subradiant state at energy ifin is 
coupled through X to the superradiant-RC subsystem (framed), where the coupling 
Vrc produces the dressed levels with energies A±. 


thus identifying the following global optimization condition; 

Xopt = V;c and = (5) 

The estimate given in ([^, obtained considering the coherent transfer of excitation 
between the subradiant and the RC states, is a very good estimate also of the global 
optimum of the transfer efficiency, as shown in hgurej^ (blue cross in left panel). 

We observe that the condition = 0 is not necessary for the probability of being 
on the RC to reach 1, but, in combination with X = Vj-c, it makes such transfer the 
fastest. On the other hand, if we have some constraint on the coupling X enforcing the 
condition X ^ Rrc, the optimal detuning is not 6E^c = 0. 

To hnd the resonant detuning producing the optimal transfer in the case X !:§> W 
we can consider R-c as a perturbation, obtaining the picture illustrated in panel (a) of 
hgure 1^ The subradiant and superradiant states couple and give rise to the dressed 
energy levels 


^±{X) 




( 6 ) 


which reduce to £±(X) = Ei^EX for A = 0 (recall that we considered both n and Ta as 
small perturbations, that can be neglected in hnding the dressed energies). The initial 
excitation is equally distributed on those levels, and we can then identify two optimal 
detuning values by the symmetric resonant tunneling conditions 


Ere = £±(^) , 


(7) 


entailing 

A /A2 

5Erc(X) = -±y—TAX (8) 

For the case A = 0 we have (5Erc(X) = EX (see dashed curves in hgure |^. Note that 
in (|^ we kept the dependence on A, since it will be relevant in the following. 
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Figure 4. Average transfer efficiency g at fluorescence time tfl = S/Ffl, computed 
starting from the subradiant state \SUB), plotted in the left panel as a function of 
the subradiant-RC detuning 5E^c and of the deterministic subradiant-superradiant 
coupling X, and in the right panel as a function of SE-^c and of the rescaled disorder 
strength W = Wjs/^. The blue cross Indicates the estimate ^ for the optimal 
transfer conditions. The dashed curves indicate the resonances determined In Q and 
the dot-dashed line marks the the resonant condition (11). The value of the parameters 
are A = 1 cm“^, V^c = 10 cm“^, n = 0.01 cm“^, and Tfl = 10“^ cm“^. 


2.1.2. The finite-gap case We will now investigate whether the optimal conditions 
given in eqnation ([^ are valid also for hnite valnes of the energy gap A. We will show 
that condition (|^ gives a good estimate for the global optimnm whenever Id-c is of the 
same order of A in the deterministic model, while disorder prodnces some modihcations, 
as discnssed in the next snbsection. 

Let ns hrst consider the situation when the energy gap A is hnite and A < Xf-c. 
We see from hgure that the global optimization condition (|^ (blue cross) is still an 
excellent estimate for the conhguration with maximal efficiency, and also the symmetric 
resonances present for X Xf-c follow the analytic prediction ([^ (see dashed curves 
in hgure 1^. Nevertheless, hgure enlightens a somewhat unexpected feature: if we 
assume now the coupling X to be constrained within the region A A < Id-c, the 
RC energy producing the maximal efficiency, identihed by a sharp resonance, is very far 
from either E[^ or e±. 

To understand such a resonance, we can now consider A as a small perturbation, 
exploiting the picture illustrated in panel (b) of hgure The superradiant and the RC 
states couple to give the dressed energy levels 

A±(^rc) = ^ (^in + A + Ere) ± (^in + A - , (0) 

which clearly depend on the RC energy Ej-c. The initial excitation is all on the subradiant 
state, since A is small, and a resonant tunneling criterion would now require Ei^ 
to match the energies X± of the superradiant-RC subsystem. Nevertheless, we have 
Ein < Ein -b A < A+ by construction, so that the resonant condition for A -C A < V(c 
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entailing 


Em — ^-{Erc) ; 


1/2 

5E,c = -f , 


(10) 


( 11 ) 


which corresponds exactly to the numerical results (see dot-dashed line in figure]^. 

If we now decrease further the ratio Id-c/we find the following remarkable result 
(see figure left upper panel): the estimate ([^, which was obtained for If-c 3> A, 
still identifies the global efficiency optimization in the deterministic model (blue cross). 
Moreover, the resonances predicted by ([8| for A S> A and by (11) for A -C A 


still correspond to the local optimization of the efficiency (see white curves in figure]^. 
As for the model in presence of disorder (right upper panel of figure |^, while the 
estimate ([^ is still within a region of significant efficiency, the optimal condition is 
modified by disorder, as a more detailed analysis shows (see lower panels of figure and 
the discussion in the next subsection). 

We finally observe that the width in energy SE,.c of the high-efficiency region at 
the optimal value of A is proportional to the coupling W, providing a measure of the 
robustness of transfer to perturbations around the optimal detuning. 


2.2. The effect of disorder 

As can be clearly seen from figuresandwhere W S> A, the results obtained from the 
deterministic model are in good agreement with those in presence of disorder. Indeed, 
we can obtain an excellent estimate for the global optimization condition by simply 
substituting the deterministic coupling A with W/y/2 in (|^ (blue crosses in figures]^ 
and|^ right panels). 

On the other side, when W ^ A (figure |^, disorder induces some modification 
of the global optimization condition. This effect can be clearly seen by comparing the 
two lower panels of figure]^ which describe a situation with Vrc/A ~ 0.5, similar to the 
one we will find in the realistic model considered below. By comparing the two lower 
panels of figure one can see that the average over disorder shifts the optimal detuning 
from 5E^c = 0 to the low-disorder resonance 5E^c = given by (0. This can 

be explained from the fact that the random coupling falls for many realizations in the 
region x < ^ A, where the resonance is for dEj-c = V^^/A and not for dE^-c = 0. As 

far as optimal disorder is concerned, even if ([^ overestimates its actual value, it still 
gives an estimate within 5% of the maximal efficiency, see figure lower right panel. 

Note that the case I^-c "C A has not been discussed since it is not relevant for the 
realistic model considered in the next section. 

Summarizing, we predict the global optimization of the average transfer efficiency 
T]{t > h/Tff) from the subradiant state of the trimer into the sink placed at the RC as 
given by 


^Vpt/v^ and ^ -E^in for U > , 


( 12 ) 
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Figure 5. Average transfer efficiency g at fluorescence time tn = h/Tn, computed 
starting from the subradiant state \SUB), plotted in the left panels as a function 
of the subradiant-RC detuning SE^c and of the deterministic subradiant-superradiant 
coupling X, and in the right panels as a function of SE^-c and of the rescaled disorder 
strength W' = lF/-\/2- The blue cross indicates the estimate ([^ for the optimal 
transfer conditions. In the upper panels the dashed curves indicate the resonances 
determined in (|^ and the dot-dashed line marks the the resonant condition (111. In the 
lower panels we zoomed on the high-efficiency regions to show how disorder affects the 
accuracy of the estimates ([^ and (12). The blue dots mark the modification (131 of the 
estimate proposed in section |2.2[ Isolines enclose the regions of efficiency lying within 
1% and 5% of the maximal efficiency. The value of the parameters are A = 20 cm“^, 
Vrc = 10 cm“^, K = 0.01 cm“^, and Ffl = 10” 


'' ^ cm 


and 

and E,, = Ei, + ^ for Kc < A • (13) 

Onr resnlts also confirm the expectation that the width in energy of the high- 
efficiency region at the optimal disorder strength is proportional to the snperradiant 
conpling Ifc- 

Another physically relevant qnestion concerns the optimal detnning at some fixed 
disorder strength determined by physiological conditions (natnral systems are nsnally 
snbject to a definite range of static disorder). In this sitnation onr resnlts indicate that, 
if the disorder is constrained to be smaller than the energy gap A and the conpling Vrc, 
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the optimal subradiant-RC detuning is not zero, but it is always given by 

1/2 

SE.C = -f , (14) 

with a signihcant efficiency present only in a narrow band around such optimal detuning. 
This result is at variance with the intuitive expectation that the best transport would 
be obtained at resonance with the initial state, 60,.^ = 0. Moreover, it can be 
easily extended to iV-sites rings and its general validity will be the subject of future 
investigations. 

For large values of the disorder, as a result of the averaging procedure, we have a 
broad resonance centered around 6Ej-c = 0 , which fades into two very broad resonances 
centered around Ej-c = e±{W/\/2) given by condition ([^, characterized by a negligible 
efficiency. 


3. Application to natural complexes 

Let us analyze a model of the light-harvesting complex I (LHI) and the RC of purple 
bacteria. Our purpose is to determine the optimal efficiency as a function of both the 
energy detuning of the RC with respect to LHI and the disorder strength. We will use 
the previous hndings concerning the trimer model to estimate such global optimum. 


3.1. LHI-RC complex of purple bacteria 


The LHI-RC complex of purple bacteria can display a variety of forms depending on 
the species (see [32] for a review). To illustrate the relevance of our analysis in such 
complexes we focus on the model presented in [3T|. It consists of 32 chromophores 
arranged to form the LHI ring with equal energies e^, plus 4 central chromophores which 
represent the RC: a special pair with energy (labeled as 33, 34) and two accessory sites 
with energy (labeled as 35,36), see £gure[^ The chromophores are modeled as two- 
level systems (sites) and, since in physiological situations the complex is irradiated by 
a dilute light, we can safely apply the single-excitation approximation. The interaction 
among sites is due to the transition dipole moments of the chromophores, except for 
the nearest neighbors where the dipole approximation is no longer valid. We then 
assume, following |3I], such interaction strength to take the values Hi,H 2 alternating 
each site, which can be determined by htting the energy spectrum of the LHI complex. 
Moreover, the two inner sites of the RC with energy e* have also a short-range coupling 
Hsp determined in a similar way. 

For positions Rg and orientations fig of the dipole moments as well as for the 
nearest neighbor interactions (Hi, H 2 , f^sp) and the site energies {er,eg,ea) we took the 


data (reported in Table Al and Table A 2 in the Appendix) from [3T] . 

The Hamiltonian of the system, written on the site-basis |s ), s = 1,..., 36, reads 


32 


34 


hfrhod = ^ |s) (s I -f ^ |s) (s 


“h Gq 


s=l 


s=33 


36 

• El 

s=35 


|s) (s 
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Figure 6. Spatial structure of the LHI-RC complex. The two red sites are the special 
pair (33, 34) (“open” sites) from which the excitation can escape. The other two sites 
in the center are part of the RC, but they are “closed”, in the sense that they contribute 
to the interaction with the other sites but the excitation can not escape from them. 
The arrows indicate the directions of the dipole moments. 


16 

+si.5:(i2*- l)( 2 s| + | 2 s)( 2 s-l|) 

15 

+ 122 ) {2s + 11 + |2s + 1) (2s I) 

s=l 

+ 122(132) ( 1 1 + | 1 ) (32 I) + ap(|33) (34 | + |34) (33 |) 

± i-gf. ± Lgj. 

where Rgr is the vector connecting the site s to r, fig is the dipole moment at the site 
s, and the constant C embeds the dipole strength common to all of the chromophores. 
The hat over the last summation indicates that we consider the interaction among all 
sites (ring and RC) but the terms associated with nearest neighbors along the ring and 
the special pair. 

The loss of excitation into the RC is modeled by two independent sinks attached 
only to the sites of the special pair (33, 34) with identical decay widths 7 rc. Hence, the 
effective non-Hermitian Hamiltonian for the system is 

34 

// = //rhod-*^ 5 ^k)(s| . (16) 

s=33 

The geometry of the system is illustrated in hgure and we want to stress the fact 
that, while a symmetry of LHI under cyclic permutations of dimers can be devised, the 
geometry of RC sites has no symmetry whatsoever. 

Even in this realistic complex it is possible to identify LHI states that are 
superradiant/subradiant with respect to the transfer of excitation towards the RC states, 
where the excitation can be trapped. To this end, we first perform two independent 



s.r 
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Figure 7. Schematic representation of the 32 energy levels of the LHI complex (left) 
and of energies and decay widths, in cm“^, of the four RC eigenstates (right) of purple 
bacteria. 


diagonalizations of the LHI and RC sectors of the total Hamiltonian. Since the LHI 
sector of H is Hermitian, it can be diagonalized by means of a real orthogonal matrix 
R'"™ (which is unitary), while the RC sector, being a non-Hermitian symmetric matrix, 
is diagonalized by a complex orthogonal matrix (which is non-unitary). By 

performing the transformation (R^^^ (g) g) R^^), the total Hamiltonian 

takes the form 


where 


/diag(£'™) a b c d \ 



a* 


0 

0 0 



b^ 

0 


0 0 




0 

0 

0 



\ i 

0 

0 

0 y 


< < ... 

E^y^, are 

the 

real) energy levels of the LHI system and 

£;RC _ ^RC 

V 

k = 

1,...,4, 


are the complex eigenvalues of the RC system, with energy and decay width Lfc. 
The vectors a, b, c, and d represent the coupling between the LHI eigenstates and the 
RC eigenstates. The intensity of such couplings ranges typically from 0 to 50 cm“^, 
with average intensity of about 6 cm“^. The energy levels of the LHI-RC system thus 
obtained are illustrated in hgure 
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By analyzing the conpling vectors a, b, c, and d, we can observe that, due to the 
peculiar symmetry structure of the system at hand, there is a great inhomogeneity in 
their intensity. Indeed, for any RC level we can identify a set of superradiant LHI 
states (strongly coupled) and a set of subradiant LHI states (weakly coupled). Such a 
distinction will be important in making the transport properties amenable to an analysis 
based on the trimer model discussed in section [2l 

The effect of the static disorder has been taken into account by adding to H the 
disorder Hamiltonian 

36 

D = ^6s\s) {s \ , (17) 

S = 1 


where the d^’s are independent Gaussian random variables centered in 0 and with 
variance W^. The disorder will produce a random coupling between subradiant and 
superradiant states with variance ~ W'^/N, where = 32 is the number of LHI sites. 
Indeed, the diagonalization of the LHI sector of (17) produces a. N x N full random 
matrix with Gaussian entries centered in 0 and with variance ~ W'^/N, whereas the 
diagonalization of the RG sector (which has 4 sites) produces a 4 x 4 full random matrix 
with Gaussian entries centered in 0 and with variance W^A. 

It should be noted that the value of the LHI site energy can be nicely determined 
by fitting the absorption spectra of the LHI complex, whereas the RG energies e* and 
Cq are reported to be lower than e^, but their precise value is not well known and they 
can be shifted to match other experimental findings such as the high values of efficiency 
(see the discussion in [3T]L 

The coupling with the electromagnetic field, which induces superabsorption and 
superflorescence of light, has been taken into account adding to the Hamiltonian in (16) 
a non-Hermitian part, following the procedure described in [TO] . 


In the following we will keep all the parameters of the model fixed, except for the RG 
energies and the strength W of the disorder, with the aim to find the global optimum. 
To achieve this goal, we will show how the main features of the transport properties of 
the LHI-RG complex are obtained by reducing it to two independent trimers. 


3.2. Reduction procedure 

In this complex we have many LHI states coupled to four RG energy levels, so that it 
may seem very difficult to optimize the transfer efficiency by choosing the RG energies 
and the disorder strength independently of the LHI state initially excited. Nevertheless, 
in physiological conditions the ground state and the first two degenerate excited states 
of the LHI ring are the most populated ones (together they gather 77% of the total 
probability at T = 288 K). Moreover, the only state of the RG we need to consider 
is the lowest energy state, since the other three RG states are separated from it by a 
significant energy gap, which makes their effective coupling to the most populated LHI 
states negligible. 
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In the doubly degenerate subspace of the hrst two excited LHI states, if we consider 
the coupling to the RC, it is possible to identify a subradiant state and a superradiant 
one, denoted by |sub) and |sr), respectively (see £gure[^. Those states are such that 

(sub I if |rc) = 02 = 0 , and (sr | if |rc) = 03 = —43.71 cni“^. 

Moreover, the LHI ground state |gs) is subradiant with respect to transfer towards |rc), 
since the coupling intensity 

I (gs I if |rc) I = fyi I = 0.26 cm”^ 

is much smaller than the average coupling, and negligible if compared to fyal = 
43.71 cm“^. 

As already noted, we view the energy value Ef-^ of |rc) as an adjustable parameter, 
which can be modihed by shifting all of the RC energies e^, of the same amount. Note 
that the hrst two excited states are optically very active [31] , concentrating most of the 
LHI huorescence, while the |gs) is optically almost inactive. On the other hand, the 
transfer from |sr) towards the RC is very fast, due to its strong coupling. It is then 
capable to compete with (and outperform) the decay due to huorescence. Moreover, 
when the coupling between |sub) and |sr) induced by disorder is strong enough, also 
the transfer from |sub) becomes significant. For this reason, in our numerical simulations 
of the LHI-RC complex we took into account huorescence in an exact way, following 
the approach presented in [10]. Such approach includes the ehect of superradiant loss 
of excitation into the electromagnetic held, thus allowing to estimate the competition 
between superradiant loss and superradiant transfer of the excitation. 

The problem of excitation transfer from an initial Gibbs populations distributed on 
the three lowest-energy levels towards the lowest RC energy can be decomposed in two 
independent problems, involving two trimers: (A) LHI ground state |gs), superradiant 
|sr), and |rc); (B) subradiant excited state |sub), superradiant |sr), and |rc). 

Concerning the trimer A, we neglect the direct coupling oi since it is much weaker 
than 03 , and we use the model Hamiltonian (|^ with the identihcations 

A = = 94 cm"^ and fyc = fyal = 43.71 cm"^. 

As for trimer B, we can use again the model Hamiltonian (|^ with A ~ 0, since 
= ^lhi ^ ^^429 cm"^, 

and fy-c = fysl- In both models the opening strength k in ([^ is given by Ti = 1.58 cm“^. 
Moreover, the variance of the random coupling between subradiant states and |sr), which 
will be important in what follows, is of the order of W‘^/N (with N = 32). 

3.3. Optimality conditions 

At room temperature the excitation is evenly shared among the states |gs) (33%), 
|sub) (22%) and |sr) (22%). To hnd the optimal conditions we will focus on the 
subradiant states, which are the only ones to display a non-monotone behavior with 
disorder. Moreover, we can use the conditions found in section for the trimer model. 
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Transfer from |gs). Given that |ai| loal, the indirect connection between |gs) and 
|rc) through |sr) is the most relevant. We then expect that, to find the optimal value 
of relative to transfer from |gs), one must consider the trimer A, which predicts 
the optimal resonant configuration at 

Efc _ ^ gLHi + « 11355 cm-', 

since Id-c/A 0.5 (see conditions (HD- Such a configuration can be produced by setting 
the energy of the special pair = Ef-^ + 1188 cm“^ = 12543 cm“^. The robustness 
with respect to detuning of such optimum has been discussed in section indicating 
that should he within an interval ild-c around Ef^^. As explained in section]^ the 
optimal disorder is given by the condition Ugs ~ W/y/N ~ Id-c = loal, which produces 
the estimate ITopt — 250 cm“^. The optimal transfer conditions for trimer A are then: 

W ~ 250 cm“^ and = E^^^ ± — (11355 ± 44) cm“^. 


Transfer from |sub). In this case (A = 0) the optimal detuning is zero (see 

± Vrc, and the optimal disorder is again given 
So that the optimal conditions for this case are 


conditions (12)), entailing Ef^ = E^^^ 


by asub ^ W/y/N ~ G, 

W ~ 250 cm“^ 


and Ef^ = Ef^ ± I/.c ~ (11429 ± 44) cm 


-1 


Global conditions. From the discussion above we see that the global optimization 
conditions are somewhat different for each of the two LHI states we analyzed, but 
the cooperatively enhanced superradiant coupling \fc is large enough to allow for a 
significant overlap of the high-efficiency regions pertaining to each initial state. This 
overlap identifies the interval 

Ef'^ ~ (11392 ± 7) cm-^, (18) 

corresponding to = (12580 ± 7) cm“^. On the other side, from the above discussion 
it is clear that the optimal disorder range is IF ~ 250 cm“^. 

This estimate of the global optimal conditions is shown as a blue cross in figure 
and it is very close to the region of maximal efficiency found numerically. Note that, 
while the detuning energy is in excellent agreement with the numerical results, the 
optimal disorder slightly overestimates the actual maximal value. This can be due to 
the fact that the efficiency of the superradiant state (not considered here) decreases with 
disorder, thus producing a lowering of the optimal disorder value. Another possible cause 
of deviation might be that for the ground state the optimal disorder can be lower than 
what we estimated, as emerges from the discussion about figure lower right panel. 
Note that the optimal disorder strength falls within the range of the natural disorder 
strengths (30 cm“^ to 300 cm“^ [23| IT2]). 

The results presented in figure also show that, for small disorder strength, the 
transfer of excitation from the LHI ring into the RC has an optimal efficiency when 
_ j^gs^ (horizontal dot-dashed line), consistent with our analysis, whereas both the 
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Figure 8. Average transfer efficiency at 20 ns for the LHI-RC complex, starting from a 
Gibbs distribution of LHI eigenstates at temperature T = 288 K, plotted as a function 
of the disorder strength W and of the energy Cs of the special pair. The horizontal 
dot-dashed line indicates the value of = 12543 cm“^ for which (resonant 

configuration). The blue cross marks the estimated optimal conditions, see text. The 
values of such that Ef-'^ = Ef™, j = 1, 2,3, (possible resonances with the first three 
LHI eigenstates) are marked with cyan bars. 

configurations with or Ef-^ = (marked with cyan bars in 

hgure would produce a lower efficiency for small W. Note that we computed the 
transfer efficiency at t = 20 ns, a time larger than the average fluorescence time. 


4. Conclusions 

We considered both paradigmatic and realistic models of quantum networks in which 
a ring of peripheral sites is coupled to some inner sites, where the excitation can be 
trapped in a reaction center. In similar networks the ring states can often be classihed 
as superradiant or subradiant, based on how they transfer the excitation into the reaction 
center. 

Subradiant states are not directly coupled to the reaction center, but static on¬ 
site disorder can effectively couple them with superradiant states. This opens an 
indirect path for the transfer of excitation from subradiant states to the reaction center 
mediated by the superradiant state. The static disorder which activates the transfer 
from subradiant states, when too strong, hinders transport so that an optimal disorder 
condition can be determined. 

We identify the building block of such kind of transport: a trimer with a subradiant 
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initial state, a target RC state, and a superradiant state coupled to the RC and, via 
static disorder, to the subradiant state. Four parameters determine the different regimes 
in which the trimer can operate: the energy distance A between the subradiant and the 
superradiant state, the intensity W of the random coupling between those states, the 
direct coupling Id-c between the superradiant state and the RC, and the detuning 5E^c 
between the subradiant initial state and the RC. 

We study how to optimize the energy transfer efficiency by varying both the disorder 
strength W and the detuning dE-^c- We determine such optimal conditions as given by: 
(i) IV ~ Kc and dErc = 0 for Kc > A; (ii) W ~ Kc and dErc = VfJA for Kc < A. 
Note that the optimal disorder is determined by the superradiant coupling V^c- We also 
discuss how to determine the optimal detuning when the disorder strength is constrained 
to be much smaller or much larger than both A and R-c- 

Notably, our general results can be applied to a realistic model for the light¬ 
harvesting complex (LHI) and reaction center (RC) of purple bacteria. A reduction of 
the whole systems to a set of independent trimers allows for a detailed explanation of the 
whole process of transfer from an initial Gibbs thermal state and justifies the very high 
efficiency of excitation transfer even at room temperature. Our analysis points towards 
a possible correlation between the structure of such a system and transfer optimization. 
Indeed, the superradiant coupling generated by the symmetry of the LHI complex is 
strong enough to allow for a simultaneous optimization of transfer from the relevant 
LHI eigenstates. Moreover, it entails an optimal disorder strength for the transfer 
from such Gibbs distribution of initial states that falls within the estimated range of 
physiological disorder. We believe that our analysis covers those quantum networks in 
which a subradiant subspace (also known in literature as trapping-free subspace [B]) is 
present. Especially to understand the enhancement in transport due to disorder, which 
is ubiquitous in such structures, it is necessary to go beyond the standard perturbative 
analysis. 

It will be important to consider, in a future work, the effect of different baths, 
represented by noise with short correlation time, which induce loss of coherence during 
the quantum evolution. Indeed, the different nature of such kind of noise, always present 
in physiological situations, and their interplay with static disorder may entail different 
estimates for the optimal noise strength, while we expect the resonant tunneling criterion 
to remain as presented in this paper. 
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Table Al. Energy and coupling parameters of the model for the LHI-RC complex. 


Description 

Symbol 

Value 

energy of the ring sites 


12911 cm-i 

energy of the special pair (33, 34) 

Ss 

12748 cm-i 

energy of the accessory pair (35, 36) 


12338 cm-i 

coupling within a dimer (2s + 1, 2s) 

ill 

806 cm“^ 

coupling across two dimers (2s, 2s + 1) 

fl2 

377 cm-i 

coupling within the special pair (33, 34) 

ilsp 

1000 cm“^ 

dipole strength of all of the chromophores 

C 

519310 A^cm-i 


Appendix A. Data for the realistic model 


The values of the parameters entering the dehnition of Hamiltonian (15) for the realistic 
model of the LHI-RC complex of the purple bacterium Rhodobacter Sphaeroides have 


been taken from |3l] and are reported in Tables Al and A2 
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Table A2. Position and transition-dipole direction of the LHI-RC sites. 


Site 
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Py 

Pz 
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